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Ultrahigh energy cosmic rays as heavy nuclei from cluster accretion shocks 
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Abstract: Large-scale accretion shocks around massive clusters of galaxies, generically expected in hier- 
archical scenarios of cosmological structure formation, are shown to be potential sources of the observed 
ultrahigh energy cosmic rays (UHECRs) by accelerating a mixture of heavy nuclei including the iron 
group elements. Current observations can be explained if the source composition at injection for the 
heavier nuclei is somewhat enhanced from simple expectations for the accreting gas. The proposed pic- 
ture should be testable by current and upcoming facilities in the near future through characteristic features 
in the UHECR spectrum, composition and anisotropy. The associated X-ray and gamma-ray signatures 
are also briefly discussed. 
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Introduction 

The origin of UHECRs with energies lO^^-lO^" 
eV and above remains one of the biggest mysteries 
in physics and astrophysics [1]. Only a few types 
of astrophysical objects appear capable of accel- 
erating UHECRs to the highest observed energies, 
such as the jets of radio-loud active galactic nuclei 
or gamma-ray bursts [2, 3]. However, no unam- 
biguous identification with any kind of source has 
been achieved so far 

In the currently favored picture of hierarchical 
structure formation in the CDM cosmology, all 
massive clusters of galaxies should be surrounded 
by strong accretion shocks, as a consequence of 
continuing infall of dark matter and baryonic gas 
[4]. Such shocks should be interesting sites of par- 
ticle acceleration, and have also been proposed as 
sources of UHECRs [5]. Here we summarize our 
recent work on this subject invoking UHECR nu- 
clei; more details can be found in Ref. [6]. 



Model 

For clusters of mass M, the rate of gas 
kinetic energy dissipation through accre- 
tion shocks can be estimated as Lace — 
9 X lO''5(Af/lOi^M0)5/3gj.gg-i [7]^ Taking the 

number density of clusters with M > lO^^Af0 
to be ris = 2 X 10~^Mpc~'^ as observed lo- 
cally within 120 Mpc [8], we can evaluate the 
total power density from cluster accretion to 
be Pace ~ iacc^s = 2 X lO^^crg s-^Mpc-^, 
which should reasonably accommodate the energy 
budget of UHECRs. 

However, estimates of the maximum energy -Emax 
for protons seem to fall short of 10^" eV by 1-2 
orders of magnitude [9, 10]. A fiducial cluster of 
A/ = 2 X IQI^A/q has shock radius R^ ~ 3.2 Mpc 
and shock velocity y^ = (4/3)(GA//i?s)^/^ ^ 
2200 km/s. The shock magnetic field is taken to 
be Bs = 1/^G, as suggested by some recent ob- 
servations [11]. The timescale for shock acceler- 
ation of particles with energy E and charge Z is 
tacc = 20k(L;)/T/2 = (20/3)(L;c/ZeB,V;2), as- 
suming the Bohm limit for the diffusion coefficient 
k{E), which can be induced by CR wave excita- 
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tion [12] and is compatible with observations of 
supernova remnant shocks [13]. To be compared 
are the energy loss timescales for photopair and 
photopion interactions with the cosmic microwave 
background (CMB) [14], the escape time from the 
acceleration region icsc ~ ^s/5k(£') [15], and the 
Hubble time tn- As is clear in Fig.l, for protons 
^max ~ 10^^-10^^ eV, confirming previous find- 
ings. 
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sider EGMF models that trace large-scale structure 
[16], as well as the case of negligible EGMF, al- 
though Galactic fields [17] are not included. A 
fraction /cr of the accretion luminosity Lace is 
converted to cosmic rays with energy distributions 
oc E-°' exp(-^/^max), and we set Ema^/Z = 
5 X 10^^ eV, a fair approximation to estimates for 
each species obtained by comparing timescales as 
in Fig.l. For the elemental composition at injec- 
tion, the He/p ratio is taken to be 0.042. All heav- 
ier elements are assumed to have the same rela- 
tive abundances at fixed energy/nucleon as that of 
Galactic CR sources at GeV energies [18], and 
scaled with respect to protons by the metallicity 
C, of the accreting gas. We take C — 0.2 as sug- 
gested by both observations and theory for the gas 
flowing in from large-scale filaments [19]. An ad- 
ditional factor A^ for the injected abundance of 
nuclei with mass number A is introduced to take 
account of possible enhancement of heavier nu- 
clei due to nonlinear modification of shock struc- 
ture by CRs [20], which may possibly be stronger 
here than for Galactic CRs due to the acceleration 
to much higher energies. 



Figure 1: Comparison of timescales at cluster 
accretion shocks for shock acceleration tacc (di- 
agnonal lines), and energy losses from interac- 
tions with background radiation fields (curves), for 
protons (thick dotted). He (thin dotted), O (thin 
solid) and Fe nuclei (thick solid). The photopair 
timescales are denoted separately for p and Fe (dot- 
dashed). Also indicated are the Hubble time tn 
(dashed) and the escape-limited E^max (circles). 

On the other hand, heavy nuclei with higher Z 
have correspondingly shorter tacc, and Fe may 
be accelerated up to lO'^" eV in the same con- 
ditions, notwithstanding energy losses by photo- 
disintegration and photopair interactions with the 
far infrared background (FIRB) and CMB (Fig.l). 
In order to explore whether nuclei from cluster 
accretion shocks can provide a viable picture of 
UHECR origin, detailed propagation calculations 
of UHE nuclei above 10^^ eV are undertaken, fol- 
lowing energy losses in the CMB and FIRB [14] 
and deflections in extragalactic magnetic fields 
(EGMF) for all particles including secondary nu- 
clei arising from photodisintegration. We con- 



Results and Discussion 

Fig. 2 shows our results for the observed spectrum 
and composition for a = 1.7 and /3 = 0.5, which 
are consistent with the current data for HiRes [21] 
and Auger [22] within about 2 sigma of combined 
cosmic and statistical variance, for both the dif- 
ferential and integral fluxes (and possibly AGASA 
[23] as well). Values of a < 2 are naturally ex- 
pected at the high energy end for nonlinear shock 
acceleration that accounts for the dynamical back 
reaction from CRs [24]. The spectral steepening at 
> 10^" eV is due both to propagation losses and 
the -Emax limit at the source. Normalization to 
the observed flux and comparison with the avail- 
able accretion power for M > 10^^ Mq fixes /cr, 
which is ~ 0.005 — 0.3 for cases with EGMF and 
~ 0.002 for the case without. Low values of fcR 
may reflect inefficient escape of CRs from the sys- 
tem, which is conceivable in view of the converg- 
ing nature of the accretion flow. CR escape may 
be mediated by diffusion in directions away from 
the filaments, or possibly by advection in outflows 
during merging events. 
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though this prediction is subject to uncertainties in 
the EGMF and Galactic fields. 
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Figure 2: Observed UHECR spectrum (top) and 
mean mass composition (bottom) versus energy E 
(1 EeV = lO^'* eV) from cluster accretion shocks 
for a ~ 1.7 and /? ~ 0.5, compared with the cur- 
rent data for HiRes (bars) and Auger (asterisks). 
The histograms are the average result over different 
model realizations for the cases with (thick) and 
without (thin) EGMF, and the thin curves outline 
the median deviations due to cosmic variance for 
the former case only. The straight line in the top 
panel denotes the injection spectrum. 

The mass composition at < 3 x 10^^ eV is pre- 
dominantly light and consistent with HiRes reports 
[25], while the rapid increase of the average mass 
at higher energies is a clear prediction of the sce- 
nario to be tested by the new generation experi- 
ments (and is in line with the latest Auger results 
on the shower elongation rate [26]). 

Despite the relative rarity of massive clusters in 
the local universe, strong deflections of the highly 
charged nuclei in EGMF allow consistency with 
the currently observed global isotropy (Fig. 3). On 
the other hand, with a sufficient number of accu- 
mulated events, clear anistropies toward a small 
number of individual sources should appear, al- 
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Figure 3: Angular power spectrum C{1) of 
UHECR arrival directions above 4 x 10^^ eV 
versus multipole I, for a realization with EGMF 
and a single, dominant cluster at D ~ 50 Mpc. 
The crosses are for 100 events with AGASA + 
SUGAR exposure and diamonds for 1000 events 
with Auger North + South exposure. Vertical bars 
indicate statistical errors. 

An aspect of this scenario that warrants further 
study is the spectral domain < 10^^ eV and the 
implications for the Galactic-extragalactic transi- 
tion region [18]. 

X-ray and Gamma-ray Signatures 

If cluster accretion shocks are indeed accelerators 
of UHE particles, we may look forward to very 
unique X-ray and gamma-ray emission that can 
serve as valuable multimessenger signals. Protons 
accelerated to 10^^-10^^ eV in cluster accretion 
shocks should efficiently channel energy into pairs 
of energy 10^^-10^^ eV through interactions with 
the CMB, which then emit synchrotron radiation 
peaking in hard X-rays and inverse Compton radi- 
ation in TeV gamma-rays. Fig. 4 displays the pre- 
dicted spectra for a Coma-like cluster, conserva- 
tively assuming that UHE proton injection contin- 
ued only for a dynamical time ~ 2 Gyr (see Ref. 
[10] for more details). The detection prospects 
are promising for Cerenkov telescopes such as 
HESS, VERITAS, CANGAROO III and MAGIC, 
and hard X-ray observatories such as Suzaku and 
the future NeXT mission. Photopair production by 
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Figure 4: Spectra of UHE proton-induced pho- 
topair emission from the accretion shock of a 
Coma-like cluster, for Bg =0.1, 0.3 and 1 /xG. 
The sensitivities for a 1 degree extended source are 
overlayed for HESS, GLAST, Suzaku XISh-HXD, 
and NeXT HXIh-SGD. 

nuclei may also be efficient and induce further in- 
teresting signals that are worth investigating. 

Combined with such complementary information 
from X-ray and gamma-rays, detailed measure- 
ments of UHECR composition and anisotropy with 
facilities such as the Pierre Auger Observatory, the 
Telescope Array, and the future Extreme Universe 
Space Observatory should provide a clear test of 
whether the largest bound structures in the universe 
are also the largest and most powerful particle ac- 
celerators. 
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